Determinant-Gravity: Cosmological implications 
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We analyze the action J d 4 x v/det \\B g^v + C R^ v \ | as a possible alternative or addition to the Einstein gravity. 
Choosing a particular form of B(R) = VR we can restore the Einstein gravity and, if B — to 2 , we obtain the 
cosmological constant term. Taking B — m + B±R and expanding the action in 1/m 2 , we obtain as a leading term 
the Einstein Lagrangian with a cosmological constant proportional to m and a series of higher order operators. In 
general case of non-vanishing B and C new cosmological solutions for the Robertson- Walker metric are obtained. 



1. Introduction. There arc numerous suggestions in 
the literature for modification of the classical Einstein 
action of general relativity: 
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SE= lS Jd i x^g(R-2A) (1) 

where R is the curvature scalar, g — — detlg^l is the 
determinant of the metric tensor, A is the cosmologi- 
cal constant and mpi is the Planck mass. Majority of 
such attempts (at least in 4-dimensional space-time) are 
based on the same structure of the action integral with 
an addition to the Einstein term of some scalar functions 
of R and/or of combinations of the Ricci, i? M „, and/or 
Riemann, R^ af 3 , tensors (/ d A x y/g C(R, R^, R^ a f}))- 
Usually, but not necessarily, one considers quadratic 
terms in the curvature proportional to R 2 , R^R^ , and 
R^vafiR^ ■ Such terms do not introduce a scale de- 
pendent parameter or, in other words, they enter into 
the action with dimensionless coefficients. Some of such 
and higher order or even non-local terms may appear as 
a result of quantum corrections, see e.g. the book [1]. 

Such a form of the Lagrangian density, i.e. a scalar 
function multiplied by the determinant of the metric ten- 
sor, is dictated by the demand of the invariance of the ac- 
tion with respect to general coordinate transformations. 
However, this is not the only way to ensure this invari- 
ance. In fact any scalar function multiplied by a deter- 
minant of a second rank tensor would also be invariant 
with respect to the choice of coordinates [2] . 

Already in 1934 Born and Infold [3] proposed a new 
action for electromagnetism (reducing to the Maxwell ac- 
tion for small amplitudes) of the form 

Sbi = J tf^detHm 2 g^+F^W (2) 

with T^y the electromagnetic field strength. 

A generalization of the above action to gravity was 
suggested by Deser and Gibbons in ref. [4], where they 
studied the physical criteria that such a theory should 
satisfy. 

Along the same line we analyze maximally symmet- 
ric spaces where only two independent tensors g^ and 



R^u are needed to specify geometrical structure of space- 
time. Within such a hypothesis we can write the general 
covariant action as: 

Sdet = J d 4 x\J det \\Q^ V \\ (3) 

where Q^ v =Bg il v + C R^y * . Below we make a simple 
choice for the coefficients: B = B(R) = m 2 +B\R where 
to, at this stage^, is a new mass parameter, while B\ 
and C are some dimensionless constant coefficients. In 
general, B and C could be arbitrary scalar functions of 
geometrical tensors, R, i? M „, R^ap, etc. 

Somewhat similar modification with C = ydet H-R^H 
was suggested by Eddington [5], as we have found out 
from ref. [6] , but the latter was considered as a pure affine 
theory depending on connection and its derivatives 
only ( a generalization of the above results in the same 
framework is given by Vollick [7]), while we prefer to 
remain in the frameworks of metric theory in the same 
spirit as ref. [4]. 

It is non-trivial task to derive equations of motion from 
action (3), but in the case of maximally symmetric space, 
considered here, this problem is greatly simplified. One 
can show that: 

^ = -B(R)V^ + B'VR^u + {g^D 2 - D^D,) 

(B'V) + (C/2) (D 2 V^ + g^D a D p V al3 - 
Dq-D^V" - D a D v Vl) = (4) 

where B' = dB/dR, is the covariant derivative and 
D 2 =g^D v . 

The second rank tensor "P M „ is defined as 



*Such an ensemble of actions is more restricted than the one 
proposed by Deser and Gibbons [4] and does not satisfy the 
condition of absence of ghosts in the linearized version [4] , [8] . 

'The Planck mass scale enters the Lagrangian only when we 
define the coupling between matter and gravity. 
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and V = g^V^. 

Remember that, according to our hypothesis, any ge- 
ometric second rank tensor, in particular, V^ v can be 
expressed through a linear combination of two basic ten- 
sors: a\ g^ v + a 2 R^, where are scalar functions of 
various invariants (R n , Tr[R ■ ... ■ R]). 

We checked explicitly that eq. (4) satisfies, as it should, 
the transversality condition, D^{l.h.s.) = 0. 

The physical consequences of the suggested above grav- 
itational action can be studied from different angles. 

In the most conservative approach we may consider 
the new action as a higher order correction to the classi- 
cal Einstein term and restrict the new interactions with 
phcnomcnological considerations. In this case we have 
to put into the r.h.s. of the equation of motion (4), 
the usual general relativity term, SttT^ / m 2 pl — G^ v , 
where T^ v is the energy-momentum tensor of matter and 
G>„ = R^u — {l/2)g fJtv R is the Einstein tensor. In the 
most radical approach, instead, we can work only with 
the new action as a fundamental dynamical source of the 
gravitational field and may even try to introduce matter 
fields under the signs of det and the square root (work in 
progress). 

In such a perspective it is important, first of all, to 
find out the proper limiting connection to the classical 
approach, eq.(l), trying to reconstruct the Einstein Gen- 
eral Relativity limit using a special form of functions B 
and C (see below). 

The first problem can be solved if we take B(R) = 
m 2 +BiR, as we mentioned above, and expand the action 
(3) in powers of 1/to 2 

S det « J d 4 x^ [m 4 + (2 Bi+ C/2)m 2 R + (6) 

(B 2 1 +B 1 C/2)R 2 + C 2 /4 [R} 2 + ...] 

where [R] 2 = R^R^ V . The first and second terms in the 
expression above correspond to the Einstein action with 
a cosmological constant (see eq. (1)) while the higher 
order terms are the corrections to be treated perturba- 
tively. It is clear that, if we want to get rid of the cos- 
mological constant, in the radical approach, we cannot 
send naively to — > because we lose the classical limit, 
hence we have to take into account canceling matter ef- 
fect (work in progress) 

We will study eq. (4) in homogeneous and isotropic 
space-time with the FRW metric: 



ds 2 = dt 2 -a 2 (t)dr 2 



(7) 



naive cancellation mechanism consist in the addition to 
eq. (3) of the term — J d A x m 4 det ||<? M „| 



As is well known it is sufficient in this case to consider 
only the time-time component of equations of motion. 

The explicit general solution of the above equation is 
quite cumbersome (and above the length of our paper) so 
we will try to find some results in physically interesting 
simple limits. 

• For example it is easy to see that a De Sitter phase 
(a(t) = e ht ) is a solution of eq.(4) of the general 



action (3) if 



rn 



(to 2 -i{AB l +C)h 2 ) =0 



(8) 



and in order to have real h solutions, the pa- 
rameters B\ and C should lie in certain bounds. 
In particular there are solutions with fixed h 2 = 
m 2 /(V3(C + AB X )) if C + 4 Si > and m ^ but 
we found that they are instable under small per- 
turbations. On the contrary, if m = 0, there are 
solutions with generic h ^ (and we checked their 
stability for any B\,C). 

• The study of possible power law solutions of cqs. (4) , 
a(t) = t n , is easier in the scale invariant case where 
m = and B — B\R. Note that the new action 
containing det \\B\ Rg^v + C R^W can be mean- 
ingful only for a definite sign of the determinant. 
In contrast to det ||<7 M „|| the latter is not guaran- 
teed from the general principles but the sign may 
be fixed on the solutions of equations of motion. 
In fig.l we show the signature of the determinant 
versus the ratio C/B\ and n = 2/(3 + 3 w). We 
see that a well defined signature in the interval 
— 1 < w < 1 is realized only for C/B\ = —4, 0: the 
second solution is, in a sense, "trivial" correspond- 
ing to the well known quadratic R 2 — Lagrangian, 
B\ j d A x^/g R 2 , while the first solutions generates 
the traceless tensor Q^ v = Bi(Rg llu — A R^ u ) whose 
equations of motion are studied below. 
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FIG. 1. Signature of det ^Rg^ + C R^\\ for a(t) = t n 
and n = 2/(3 + 3w). The negative values are shown as black 
areas, while the positive ones are whites. 
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For the power law solutions, a(t) — t n , we find the 
constraint 
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3 + 3(2 + a;)(6 + .T) 



(9) 



where n = 2/(3 + 3 u>), x = C/B\. 

This curve is showed in fig. 2. In order to con- 
trol stability of the solutions we have found that 
the first curve to the left corresponds to unsta- 
ble solutions. For the central one there is a sta- 
bility region for < w < — 1 corresponding to 
— 2V3 < x < 2(— 3 + \/3), and for the third curve 
to the right, stability is ensured if— 1/3 < w < 
and x > 2v / 3- Thus, there is no stable solution for 
w > and the matter dominated regime (w = 0) 
is the stability border line. 

For example the solution a(t) = t 1 / 2 corresponding 
to w = 1/3 (radiation dominated universe) are ob- 
tained for x = —4 and for x = and both results 
are unstable. 




FIG. 2. Plot of eq.(9) which gives the values of w as func- 
tion of C/Bi indicating when the power law a(t) = t 2 /( 3 + 3tu ) 
ansatzs is a solution of the scale invariant Lagrangian with 
m = 0. 

• The minimal scale invariant case with m = B\ = 
Sdct -> C 2 J d 4 x^Jdet\\R^\\ (10) 

gives a tractable equation of motion: 
X 2 X 2 + 6 X + 3 




FIG. 3. Solutions of eq.(ll) in the case Hi > 0, 
where H(t) is evolving for different initial conditions 
H'JHi = 1.5, 1, 0.3, 0, -0.1, -0.5, -1 and H { = 1 . 

At expansion regime, when H > the solutions 
tend to the De Sitter one with H — ► const both for 
negative and positive initial values of H (see fig. 3), 
while at contraction, when H < 0, the solutions 
may tend to singularity, H — > — oo if Hi < or to 
H = if the initial value of Hi is positive (see fig. 4). 
However, after reaching zero value. H docs not 
change sign but becomes negative again and also 
reaches singularity. Thus, in this approximation 
singularity is not avoided, though the contraction 
may stop for a while. 
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FIG. 4. Solutions of eq.(ll) in the case Hi < 0, where H(t) 
is evolving for different initial conditions H^/Hi = 0.3, 0.5 , 3 
and Hi = -1 . 

• The scale invariant case with m — and C — — 4 Si 
also deserves some comments. The action in this 
case is: 



= (11) 



d xJ det\\Rg^„ - iR^\ 



(12) 



where H = a/a is the Hubble parameter and x = 
H/H 2 . 



The determinant of Q^ v is sign definite and the ratio 
dctQ^u/ dctg^ — (2H) 4 . The equation of motion 
is very simple: 
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2HH + 6H 2 H - H 2 = (13) 

and can be integrated analytically in quadratures 
but numerical solution is taster and simpler. At ex- 
pansion regime, when initially Hi > 0, the solution 
tends to a Dc Sitter one with H — ► const. On con- 
traction phase, when Hi < 0, we find that either H 
reaches negative and infinitely large value in finite 
time, if Hi < 0, or H tends to zero, again in finite 
time if Hi > 0. 

Surprisingly equation (13) is the same that can 
be derived from the simple quadratic Lagrangian 
J d 4 Xy/g R 2 [11] or, in our approach, from action 
(3) with m = and C = 0. 

It is interesting to study eq.(13) in presence of the 
Einstein term and matter, that means we have to 
add the term k(T q — Goo) to the right hand side, 
where k is a constant, Too = P is the energy den- 
sity of matter, and Goo = 3fi 2 . Here we took the 
system of units with m 2 Pl /&Tt = 1. 

In this case the behavior of the solution at expan- 
sion regime becomes different. If n > the solu- 
tion H(t) has an oscillating behavior. For vacuum 
energy, when p = const, the oscillations proceed 
around a constant average value (see fig. 5), while 
for normal matter with p ~ a - ™, where n = 3 for 
non-relativistic matter and n = 4 for relativistic 
matter, the oscillations proceed around gradually 
decreasing average values (see fig. 5). The accelera- 
tion parameter, (a /a) = H + H 2 also oscillates and 
may be both negative and positive, depending upon 
time moment. The line H = is a separation line 
of equation (13) and solutions can smoothly touch 
it but never cross. 
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FIG. 5. Solutions of eq.(13) in presence of matter and 
cosmological constant (p = pmatter + pA — eT iz + 0.01) 
(z = log[o(t)]) in the case Hi = 1, H[ = 0, k = 10. We plot 
the value of logoff] and the descending straight line represents 
log[p matter ] = — 3z while the horizontal one is the attraction 
point 3H 2 = pa 



As we have already mentioned, the determinant of 
is sign definite only for C/B\ = 0,-4. In other cases 
it may change sign in the course of evolution and the 
solution would encounter quite strange singularity. To 
avoid that we may try the action with the absolute value, 
I d et 6 L | • It preserves invariance with respect to coor- 
dinate transformation, since the latter does not change 
the sign of determinant of a second rank tensor. Equa- 
tions of motion in this case acquire an additional factor 
containing the sign function, e[z] = z/\z\. In particular 
the only change in eq. (4) is the change of V^ v by the 
absolute value under the sign of the square root. 

As a result the singularity, at the point where the de- 
terminant vanishes, disappears and the equations of mo- 
tions can be solved for all values of time. We have found 
some solutions of time-time component of the equations 
of motion for a particular case of Q^ v coinciding with the 
Einstein tensor = C(ii jUI/ — g tlu R/2). 

In this particular case, the equation for time-time com- 
ponent has a very simple form: 

H = -3HH + kJ \2H + 3ff 2 | e[2H + 3H 2 } 3I **~ P (14) 
v H \H | 

where k is a constant which may be both positive and 
negative, e[z] is the sign function as has been already 
mentioned above. This sign function appears as a result 
of taking the absolute value of the determinant of Q^u, 
while | if | appears because this determinant is propor- 
tional toVW^ \H\. 

We have solved this equation for different special cases. 
At expansion regime, when H > 0, the solution without 
matter tends to H = 2/3t, i.e. to the matter dominated 
(MD) regime but without matter (see fig. 6)). 
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FIG. 6. Solutions of eq.(14) in absence of matter (p — 0) 
and with the initial conditions H t — ±1, H[ = 0, k = — 1. 
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FIG. 7. Solutions of eq.(14) in presence of the cosmo- 
logical constant, p vac = 3, and with the initial conditions 
Hi = 1, H[ = 1 for different values of n — —0.3, —0.1, —10. 
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FIG. 8. Solutions of eq.(14) in presence of matter den- 
sity (p = 0.1/a(i) 3 , a(0) = 1) with the initial conditions 
Hi — —1, H[ — and different values of k — —0.1, —1, —10. 



The same behavior is found in presence of relativis- 
tic and non-relativistic matter. If vacuum energy is 
non-zero and initially sub-dominant then the solution 
starts from MD regime and then turns into De Sitter 
one with H — ► const. An interesting feature of this 
solution is that before reaching an asymptotic constant 
value the solution oscillates around it with the frequency 
uo 2 w nm 2 pl /8ir (see fig. 7). Such fast oscillations would 
give rise to gravitational particle production with ener- 
gies near the Planck energy. This process might be effec- 
tive at the early inflationary stage and probably even now 
when H 2 m 2 Pl dropped down to the present day value of 
the vacuum energy. It is tempting to identify these grav- 
itationally produced particles with the sources of ultra- 
high energy cosmic rays. 

For the case of initially negative H (contraction 
regime) the solution in the absence of matter and vac- 
uum energies demonstrates very interesting feature that 
the Hubble parameter grows up to zero, becomes pos- 
itive and after reaching some maximum value turns to 
MD regime. Thus in the "empty" universe cosmological 
solutions always turn into expanding ones with H = 2/3t 
even if started from contraction. In this simple case cos- 
mology is non-singular and contraction turns into expan- 
sion (see fig. 6)). 

More interesting and more realistic situation is non- 
vanishing energy density of matter and/or vacuum. In 
this case the line H = cannot be crossed, at least as 
we were able to check on some simple examples (see fig. 
8). We have not found any realistic case of non-singular 
contraction but cannot exclude that such may exist. 



This approach can be also implemented in multidimen- 
sional case [9] . It is worth to note that a change in num- 
ber of dimensions does not request a change in the di- 
mensionality of the coefficients in the action. 

As any other higher derivative modifications of grav- 
ity, the theory introduced here possibly encounters all 
the problems which such theories are known to have in 
quantum case [10], [4]. Still more detailed investigation 
seems necessary. 
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